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£N| . Abstract. Under natural energy and decay assumptions, we derive a priori estimates for solutions 
of a Schrodinger-Newton type of equation with critical exponent. On one hand, such an equation 

^ \ generalizes the traditional Schrodinger-Newton and Choquard equations, while, on the other hand, 

■ it is naturally related to problems involving scalar curvature and conformal deformation of metrics. 

m 

p""! " 1 1. Introduction. 

\ We shall study the behavior of positive solutions to the equation 

Cj / 1 2n \ n+2 

-(-> ' A.u + ( -: — 7j * u n ~ 2 )u n - 2 — Vu = (1.1) 

& \ V |x| £ / 

in M. n , n > 3, where A is the Euclidean Laplacian, * means convolution, I is a real number, and V 
is a smooth real valued function. Equation (1.1) will be referred to as critical Schrodinger-Newton 
equation. We are concerned with a priori estimates, i.e., bounds on u and its derivatives, which 
are necessary conditions for any (positive) solution of (1.1). Needless to say, not only are a priori 

■ estimates one of the main tools towards an existence theory for a given PDE, but also they reveal 
t^- [ important features about the behavior of solutions. The form of the bounds one generally seeks to 

establish depends, of course, on specific characteristics of the equation, which in the case of (1.1), 



> 

o 



o 



CN . will be captured by suitable hypotheses on £, V and asymptotic conditions for u. In order to state 



natural assumptions for the critical Schrodinger-Newton equation as well as to highlight why one 
would consider (1.1) in the first place, we first turn our attention to some related problems. 
Recall that the Schrodinger-Newton equation is given by 

ih !& = _ti!_ A9 _ Gm 2f _L *N/| 2 U, x G M 3 , (1.2) 
5_i . at 2m V |x| / 

where \l/ = ^>(t,x) is a function onlx R, and m, h and G are constants. Physically, \& is the 
wave-function of a self-gravitating quantum system of mass m with gravitational interaction given 
by Newton's law of gravity; h and G are Planck's and Newton's constants, respectively. Equation 
(1.1) is obtained by considering the Schrodinger equation 

<9* fi 2 

ih— = A^ + m[/^ (1.3) 

dt 2m y 1 

with a Newtonian gravitational potential U that is sourced by a distribution of mass given by ^ 

itself, 

AU = 4vrGm|*| 2 . (1.4) 

In other words, the mass distribution is given in probabilistic terms, with its probability amplitude 
evolving according to the Schrodinger equation, as is usual in quantum mechanical systems. Writing 
U in terms of the right hand side of (1.4) and the fundamental solution of the Laplacian, and using 
the resulting expression into (1.3), formally produces (1.2). 
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The Schrodinger-Newton equation was first introduced by Ruffini and Bonazzola in their study of 
equilibrium of self-gravitating bosons and spin-half fermions [59] and gained notoriety with Penrose's 
ideas about the role of gravity in the collapse of the wave function [54, 55]. More recently, it was 
used in discussions of semi-classical quantum gravity [14, 35, 58]. 

Notice that the power of the convoluted term ^ becomes n — 2 in higher dimensions. Parallel to 
this situation, the following generalization of (1.2) has been considered, 

i^ = A^ + p(^H^M p - 2 , ^K", (1.5) 

where p > 2 and £ € (0, n). For the remainder of the paper, as in (1.5), dimensional constants such as 
h and G are set to one. Equation (1.5) is used in certain approximating regimes of the Hartree-Fock 
theory for a one component plasma; see e.g. [42, 43]. In three spatial dimensions, with I = 1 and 
p = 2, equation (1.5) has been extensively studied, see [7, 15, 32, 34, 44, 45, 57, 63] and references 
therein. 

As in many situations in Physics, one is particularly interested in wave-front-like solutions of the 
form (p(t,x) = e lult u(x), which, upon plugging into (1.5), leads to 

Au + uju + p(^-^y * \u\ p ^ju\u\ p ~ 2 = 0. 

Considering the more general situation where u i— > — V = —V{x) and dropping the factor p in front 
of the convolution, we find 

Au - Vu + (-j-p * |u| p ) u\u\ p ~ 2 = 0, (1.6) 

which is referred to as generalized non-linear Choquard equation. A detailed study of (1.6), including 

existence results, has been recently carried out by Ma and Zhao [48]; Cingolani, Clapp and Secchi 

[19]; Clapp and Salazar [20]; and Moroz and van Schaftingen [52]. These works deal with the case 

where the exponent p is sub-critical, i.e., p < while the case p = is called critical. Criticality 

is here understood in the usual sense of the Sobolev embedding theorems. We recall that, roughly 

speaking, equations with sub-critical non-linearity are suited for treatment via calculus of variation 

techniques (see e.g. [5]), provided that the equation can be derived from an action principle — which 

is the case for (1.6), see the above references. 

Equations with critical exponent appear in several situations in Physics and Mathematics (see e.g. 

[5, 18, 39] and references therein). One important case is the Yamabe equation 

n — 2 n+2 

A n u ~ 77 7^ R ci u + = 0, u > 0, (1.7) 

4(n — 1) 

where A g and R g are, respectively, the Laplace-Beltrami operator and the scalar curvature of a given 
metric g, and K is a constant. Equation (1.7) figures in the famous Yamabe problem [65]. We recall 
that this corresponds to finding a constant scalar curvature metric in the conformal class of a given 
closed 1 n dimensional (n > 3) Riemannian manifold. The complete solution of the Yamabe problem 
through the works of Yamabe [65], Trudinger [64], Aubin [6] and Schoen [61] was probably the first 
instance of a satisfactory existence theory for equations with critical non-linearity (see [41] for a 
complete overview). The analogous equation for the Euclidean metric was studied in great detail by 
Caffarelli, Gidas and Spruck [13]. 

Interest in the Yamabe problem has not faded with its resolution. On the contrary, the discovery 
of Pollack, that it is possible to find an arbitrary large number of solutions to the Yamabe equation 
on manifolds with positive Yamabe invariant, has led to an intensive investigation of the properties 



The Yamabe problem for manifolds with boundary was studied in [2, 3, 4, 10, 16, 21, 25, 26, 27, 30, 31, 36, 37, 49, 50]. 



CRITICAL SCHRODINGER-NEWTON 



3 



of the space <1> of solutions to (1.7) — see [41] for a definition of the Yamabe invariant and [56] 
for a precise statement of Pollack's result. A quite satisfactory account of the topology of $ was 
given through the combined works of Khuri, Marques and Schoen [40]; Brendle [9]; and Brendle 
and Marques [11] (see also [22, 46, 47, 51, 60, 62] for earlier results). These results imply that <J? is 
compact in the C 2,a topology for n < 24 and non-compact otherwise 2 . Such results were extended 
to manifolds with boundary in [21]. 

From an analytic perspective, the richness surrounding equation (1.7), including the surprising 
cut-off in dimension n = 24, is a direct consequence of the critical exponent. It should be expected, 
therefore, that allowing p = m (1-5) will lead to many interesting new phenomena, adding to 
the already sophisticated nature of the generalized non-linear Choquard equation. A contribution 
in this direction is the goal of the present work. In order not to lose sight of the relation between 
what has been just described and our objectives in the rest of the manuscript, notice that from the 
point of view of the theory of partial differential equations, the aforementioned compactness of <I> 
corresponds to a priori bounds for solutions of (1.7). The reader should also notice the similarities 
between (1.1) and (1.7), specially if we are given a metric g in M. n , with A g replacing A and V being 
the scalar curvature. 

We shall present a priori estimates for positive solutions to (1.1). Such estimates constitute the 
first step towards an existence theory for this equation. They also provide insight on the structure 
of the space of solutions to (1.1), at least for those solutions satisfying some additional requirements. 
We also give an account of the profile of blowing-up solutions. We point out that since, to the best of 
our knowledge, equation (1.1) has not been considered before in the literature, we shall not attempt 
to derive very general results; rather, our focus will be on conditions that allow, on one hand, a good 
grasp on the behavior of u without, on the other hand, rendering the problem uninteresting. We 
also stress that our methods may shed new light in the study of equations (1.2) and (1.5), in that we 
investigate the pointwise behavior of solutions as opposed to the 1? techniques previously employed 
to deal with these equations. 



Notation 2.1. From now on, u will denote a positive solution of (1.1). 

The first thing we investigate is the range of I values which will be allowed. For the integral 



to be finite near the origin without expecting any vanishing of u in the neighborhood of zero, we must 
have I £ (0, n). Next, we ask what kind of asymptotic behavior should be required for u. Experience 
with equations with critical exponent [13, 41] suggests that we should adopt 



Then (2.2) and £ € (0,n) guarantee that the integral (2.1) is finite. 

Definition 2.2. Given real numbers g > and L > 0, we say that u has (g, L)-decay if it satisfies 

u(x) < L\x\ 2 ~ n , for |x| > g. 
Denote by C Qj l the set of solutions u with (g, L)-decay. 

2 This under the assumption that the Yamabe invariant is positive, and the underlying manifold is not conformally 
equivalent to the round sphere. The cases of negative and zero Yamabe invariant are trivial. The geometric reasons for 
singling out the sphere, and the relation between the compactness of $ and the geometry of the manifold, are discussed 



2. Setting and statement of the results 




(2.1) 



u = 0(\x\ 2 n ) as |x| -t oo. 



(2.2) 



in [12, 53]. 
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We shall also need some energy conditions. In order to motivate them, multiply (1.1) by u, 
integrate by parts, and assume that all the integrals are finite. Then 

f n o\ f f 1 2ra 2n 

/ (\Vu\ 2 + Vu 2 )= / / —- I u^-2(x - y)u^--i{x)dydx. (2.3) 
jR n Jm. n Ju n \y\ 

The left-hand side is just the energy associated with the linear operator A — V. If we had a constant 
i 2n 

rather than -rp * u™- 2 , the above expression would produce the analogue of the Yamabe quotient 
for our equation. This motivates the following. 

1 2n 

Definition 2.3. We call the convolution q u (x) := t4j * u n ~ 2 the quotient of u (which is always 

\x\ 

non- negative) . For a given real number K > 0, denote by Qk the set of solutions u whose quotient 
is less than or equal to K. More precisely 



Qk ■-- 



q u (x) < K for all i£R"}. 



Remark 2.4. Since q u is related to the energy on the left-hand side of (2.3), u G Qk can be thought 
of as an energy-type of condition. This should not be confused, however, with the more physically 
appealing notion of energy for (1.5) used in [19, 20, 48, 52]. 

In order to motivate the extra hypotheses that will be needed, we have to say a few words about the 
general situation that will be investigated. We shall employ blow-up techniques to study sequences 
{ui} of C 2,a solutions to (1.1), and we are primarily concerned with the constraints that the equation 
imposes on blow-up up sequences. Understanding how a sequence {ut} can be unbounded in, say, the 
C°-norm, is important not only because such families of solutions are obstructions to the application 
of standard compactness arguments, but also because this type of behavior is expected for many 
critical equations (see [1, 8, 17, 23, 24, 28, 29, 38] and references therein). 

Whenever blow-up occurs, i.e., || Ui ||c°(R n ) — ^ oo as z — > oo, condition (2.2) restricts the blow-up to 
within a compact set, in which case, we can assume Ui to diverge along a sequence of points xi — > x. 
An analysis of the sequence {u^} is carried out by rescaling the solutions and the coordinates, leading 
to an appropriate blow-up model for equation (1.1). In this situation, one expects that the blow up 
of Ui, together with {u{\ C Qk, implies that the rescaled q Ui 's are very close to a constant in the 
neighborhood of x. But in order to avoid substantial extra work that would distract us from the main 
goals of the paper, we shall simply assume that q Ui has this desired property. Moreover, although 
our analysis will be local in nature, to avoid the introduction of further cumbersome hypotheses, 
such an assumption will be taken to hold on a big compact set. The precise behavior of q u has to be 
ultimately determined by a more refined analysis of the solutions to equation (1.1), what is beyond 
the scope of this paper. A technical condition on the rate at which q Ui approaches a constant will 
also be assumed, although probably this can be relaxed. With this in mind we now state our results, 
whose essence is that control over the convolution q u yields uniform control over the solutions. From 
elliptic theory, one would expect that the required control on q u should be in the C°' a topology, and 
that turns out to be in fact the case. 

Theorem 2.5. Fix positive numbers g, L and K. Let m be a sequence of C 2,a positive solutions to 
(1.1), < a < 1, satisfying 

{ Ui }Zi C C e , L n Q K , 
and suppose that there exists a constant Q such that 

q Ui ^Q>0 in C°' Q (£ r (0)) as i -> oo, 
for some r > g. Suppose further that n > 6. 
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If II u i llc°(M n ) — ^ 00 an d II Qui — Q llc*°' Q (B r (o)) II u i llcopst™) - ^ as ^ — °°; i/ien, up to a subsequence, 
the following holds. There exist x G B s (0), a sequence Xi — )• x and a positive number a such that 

II n 'i llc°(K") = Ui(xi) 

and 

|| (^(xj)) -1 ^ - (ui(xi))~ l Zi \\ c o(B tT {x))^ as i -> 00, 

waere 

-4- Q 



Z;(x) := (Ui(xi)) 1 f(ltf(xi)) "- 2 7 rr \x-Xi\ 2 ) 

\ n(n — 2) / 



4 

re— 2 



n(n 

Furthermore, the following estimate holds 

|| (^(Xi))"^- (UifXj))" 1 ^ ||cO(B CT (x))< C(Ui(xi))" 

w/iere C = C(L, g,K,Q,r,n,a, \\ V ||co.<*(.B r (o)))- 

Remark 2.6. The restriction to n > 6 is used to obtain the most direct proof without consider- 
ing variations of (2.2). We believe that this can be removed by a more careful application of the 
techniques here presented. 

Theorem 2.5 identifies a blow-up model for equation (1.1). In other words, it states that under 
suitable energy and decay conditions, and up to a subsequence, any family of solutions that blows up 
is approximated, after rescaling and near the blow up point x, by the radially symmetric functions 
Zj. The following corollary says that if we also rescale the coordinates, then the C° convergence of 
theorem 2.5 is improved to C 2 convergence. 

Theorem 2.7. Assume the same hypotheses and notation of theorem 2.5. If \\ Ui ||c7°(]R n ) — * 00 an d 
II Qui — Q \\c°- a (B r (o)) II u i \\c~°(R n )~ * as i ~ * °°> letting Xi and x be as in the conclusion of theorem 
2.5, the following holds. Define 

_-. 2 
Vi(y) := (ui(xi)) Ui(xi + («i(x»)) a -«3/) 

and 

— 1 2 
Zi(y) := (ui(xi)) Zi(xi + (u^Xj)) 2 -»?/). 



2 — n 



Then Z^y) = (1 + n{ ^_ 2 ) \v\) 2 = Z(y) for every i, and 



v% - Z ||c 2 (B A .(o))^ as i 00 ) 



where Xi = (ui(xi)) n ~ 2 rj, with r/ a small positive number. Furthermore, the following estimate holds 

4 

\\Vi~ Z \\ C l(B Xi (0)) < C{Ui{xi)) »-2, 

waere C = C(L,g,K,Q,r,n,a, \\ V ||c°.«(B r (o)))- 

The function Z in theorem 2.7 is well known in conformal geometry. The metric on the sphere 
written in a coordinate chart via stereographic projection is conformal to the Euclidean metric, with 
the conformal factor being a multiple of Z. Z has also been extensively used in the study of the 
Yamabe problem. 

One naturally wonders about the boundedness of families of solutions to (1.1) as well as the 
possibility that sequences Ui > degenerate in the limit, i.e., points where limj_ 5 . 00 Ui(x) = 0. The 
next theorem gives some sufficient conditions for bounds from above and below on u. 
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Theorem 2.8. Assume the same hypothesis of theorem 2.5, suppose that V ^ and that it does not 
change sign in B r (0). A sufficient condition for the existence of a constant C = C{L, g,K,Q, r, n,a, \\ 
V \\c°' a (B r (o))) such that the inequalities 

II u llc°(R«) + II u \\c 2 ' a (B e (0)) < C, 

and 

1 

C 

hold for any u G C q .l n Qk, is that the inequality of theorem 2.7 be true in the C°-norm with 

4 4 g 

(ui(xi)) n ~ 2 replaced by (ui(xi)) n ~ 2 , for some 5 > 0. 

3. Proof of the theorems 

The hypotheses and notation of theorem 2.5 will be assumed throughout this section. Let {ui} be 
such that || Ui \\c°(R n ) — ^ oo as i — s> oo. Since u has (g, L)-decay, we can assume that for large i 



u llc°(B e (0)) > 



II Ui ||cO(R«) = II u i llc* (Bp(0)) • 

Letting Xi be such that || Ui \\ c o^ B (o) ) = u ( x i)> an d passing to a subsequence we can assume that 

Xi — > x for some x 6 B S (Q), and taking i sufficiently large if necessary we can also suppose that Xi 
and x are interior points. Define a sequence of real numbers {ei}^ zl by 

2-n 

e ■ 2 :=Ui(xi). 

Notice that e% — > when i — > oo. £% measures the rate at which m blows- up. For each fixed i, 
consider the change of coordinates 

y = e~ l {x - Xi), 

and define the rescaled functions 

Vi(y) := e~Ui(xi + e^y). 
Then is a local maximum for Vi with Vi(0) = 1 and 

< Vi < 1. (3.1) 

A direct computation shows that Vi satisfies 

Auj + q iVz n - 2 ~ efVv, = 0, (3.2) 

where qi{y) = q Ui (xi + £iy) and V(y) = V(xi + £iy). Let Zi and Z be as in theorems 2.5 and 2.7. 
In the sequel, we shall evoke several standard estimates of elliptic theory. A full account of these 
results can be found in [33]. 

Lemma 3.1. With the above definitions, up to a subsequence, it holds that Vi — >■ Z in Cf oc (M. n ). 

Proof. Fix R > and d > 0, and set R' = R + d, R" = R + 2d. For any p > 1, we have by IP 
estimates that 

II ^ \\w 2 >p(B r ,{0))< C( || Vi \\lp(B r „(0)) + II QiV™~ 2 \\lp(B r „(0)) +£i II V'Ui 11^(5^,(0)) ), 

where C = C(n,p, R' , R") and W 2,p is the usual Sobolev space of functions with 2 weak derivatives 
in LP . From (3.1) and our hypothesis it follows that 

II Vi \\w 2 >p(B R ,{0)) < C(! + ^ 2 ) < 2 C> 
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for large i, where C = C(n,p, R' , R" , K, \\ V \\c°(B Rll (o)))- Choosing p such that 2 > ^, we then 
obtain by the Sobolev embedding theorem that Vi is bounded in C l,a ' (B R i(0)), for some < a' < 1, 
therefore there exists a subsequence, still denoted V{, which converges in C 1,a (i?/j'(0)), < a < a', 
to a limit Voq. 

Next, evoke Schauder estimates to obtain 

II v i \\c 2 < a {B R {0))< C{ || Vi \\c°(B Rl (0)) + II Qi v r~ 2 \\c°><*(B R ,(0)) + £ i II Vv i llc*°> Q (B Rl (0)) ): 

where C = C(n,a, R, R'). Combining this inequality, the interpolation inequality, the previous 
bound on the C 1,Q -norm of and the hypotheses of theorem 2.5, we conclude that 



Vi \\c^(B R (o)) 



< a 



with C = C(n, a,p, R, d, K, Q, \\ V \\c°(B R „(o)))- As a consequence, up to a subsequence and decreas- 
ing a if necessary, the above C 1,a convergence is in fact C 2,a convergence, and we can pass to the 
limit in (3.2) to conclude that satisfies 



n+2 
. n-2 



Av^ + QvSo 2 =0 in B R (0). 

Take now a sequence Rj — > oo. Using the above argument on each B Rj (0) along with a standard 
diagonal subsequence construction produces a subsequence v i that converges to a limit Voo satisfying 

n+2 

Avoo + QvSo- 2 = in R n , (3.3) 

with the convergence being in C 2,a on each fixed Br (0). Therefore v — > Voq in Cf^(M. n ). Solutions 
to (3.3) have been studied by Caffarelli, Gidas and Spruck in [13]. From their results and the fact 
^oo(O) = 1, we obtain 



2-n 



Q 

□ 



«oo(y) = (i + - r ^- 5T |y|) 2 =Z{y) 

n(n — 2) / 



Proposition 3.2. There exists a constants C > 0, independent of i, such that 

\vi - Z\(y) < Ce 2 for every \y\ < er 1 , 
possibly after passing to a subsequence. 
Proof. Set 

Ai= max \vi - Z\ = \vi - Z\(yi), 
\y\<e^ 

where yi is defined by this relation as a point where the maximum of \vi — Z\ is achieved for \y\ < s^~ l . 

Suppose first that, up to a subsequence, \yi\ < If the result is not true, then there exists a 

subsequence such that 

e 2 

-j — > as i -> oo. (3.4) 
Define Wi(y) = A~ l (vi — Z)(y). Then it satisfies 

A Wi + oiitfj = -i- [ Vvi + ^2^^"~ 2 , (3.5) 



-4; 



-7 
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where 

n+2 



n + 2 



<H = Q- ? • 

By Taylor's theorem, 

™+2 n + 2 n + 2 4 6-n „ 

v;~ - Z^ = ^—Z^(vi - Z) + 0(Z^\ Vi - Z\ 2 ). (3.6) 

From our hypotheses, expansion (3.6), the (g, L)-decay of u, lemma 3.1, (3.4) and equation (3.5), 
it follows by an argument similar to the proof of lemma 3.1, that Wi is bounded in C 2 oc (IR n ) and 

C 

- (1 + 13/1)4 ' (3 - 7) 

for some constant C > independent of i. Passing to a subsequence, we see that Wi — > if^oo in 
Cf oc (M. n ), and that Woo satisfies 

11+2 4 

A^oo + -QZ— Woo = 0. (3.8) 

n — 2 

On the other hand, using the Green's function for the Laplacian with Dirichlet boundary condition, 
the representation formula and (3.7) show that for |y| < 

\^(v)\<T^ri + C £ f, (3.9) 
1 "r \y\ *H 

for some constant C > independent of i. In particular, Wqq has the property that 

lim Woo (y) = 0. (3.10) 

\v\-*x> 

Up to a harmless rescaling, solutions to (3.8) with the property (3.10) have also been studied by 
Caffarelli, Gidas and Spruck in [13]. They have the property that i^oo = if ^oo(O) = = |VtOoo(0)|. 
Recalling that Vi(0) = 1 and that is a local maximum of Vi, and noticing that Z(0) = 1, VZ(0) = 0, 
we see that Wi(0) = = |Vu?j(0)|. Therefore, Woo(0) = = |Vwoo(0)| holds, and hence Woo vanishes 
identically. Since Wi(yi) = 1 by construction, we must have \yi\ — > oo, but this contradicts (3.9) 
because of (3.4) and Wi{yi) = 1. 

It remains to prove the proposition in the case when > Notice that from the (g, L)-decay 
of u, we obtain that Vi(y) < C\y\~^~ for \y\ > geT 1 , while Z obeys the estimate Z(y) < C(l + |y|) 2_n . 
From these we get \vi — Z\(yi) < C\yi\~^~ < Cef if \yt\ > +p. □ 
Proof of theorem 2. 5: Notice that 

2— n 

Zi(x) =e i 2 Z{e7 l x). 

Writing the estimate of proposition 3.2 in x-coordinates and recalling the definition of e%, we obtain 

|| (u^Xi^u - (uiixi))' 1 Zi || c o( Bl(a ..))< Cej. 
Since xi -+ x, choosing a > small and i large, we obtain the result. □ 



Proof of theorem 2.7: For the C°-norm this is simply the estimate of proposition 3.2 written in 
x-coordinates. From the C° bound, we obtain the C 2 bound by standard elliptic estimates, after 
considering a smaller ball via the introduction of rj. □ 



CRITICAL SCHRODINGER-NEWTON 



Proof of theorem 2.8: We shall first show that || u\ ||co(H n ) < C for a constant independent of i. If 
this is not the case, then we can find a sequence m that blows up as in the assumptions of theorem 
2.5. In light of elliptic theory and shrinking r\ if necessary, we obtain that the hypothesis on the C° 
decay of Vi — Z then yields the estimate 



z \\c*(b x .(o)) < Cel +S , 



which then gives 



n— 2 



V m ( Ui - zt) Hew**)) < Ce 2 + S ~ m , m = 0, 1,2. (3.11) 



n-2 n-2 / n+2 n+2 . n-2 , \ n + 2 



II Y7«V„,. , << m 

Notice that 

e^A( Ui - zi) + Qe^~ (<" 2 ~ <" 2 J + [q Ui - Q)u-- 2 = e^Vm. (3.12) 
From (3.11), and the hypotheses on q Ui we conclude that 

n — 2 / n + 2 n + 2 s n-2 / \ n + 2 

2 Q(< 1 - 2 - zr 2 )(^) + e t 2 ^ - QJ (xi)ur 2 (x<) as i oo. (3.13) 

(3.11) also gives 

n — 2 

£j 2 A(tij - Zi){xi) ^0 as z -> oo. (3-14) 
But (3.13) and (3.14) yield a contradiction with (3.12) since 

n-2 

e i 2 Ui(xi) = 1 for every i, 

and V is bounded away from zero in B g (0). 

This establishes a bound on the C°-norm of U{. Restricting the problem to B r (0) and using 
standard elliptic estimates, produces a bound on the C -norm within B e (0). With the C 2 bound at 
hand, we write 

n+2 4 

=<- 2 Ui, 

4 

and consider u™~ 2 as a given coefficient. In this case, we can treat the equation as a linear equation 
(for Ui) for which the Harnack inequality can be applied, producing the desired bound from below 
on Uj. □ 
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